Parametric amplification of metric fluctuations through a bouncing phase 
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We clarify tlie properties of the beliavior of classical cosmological perturbations when the Universe 
experiences a bounce. This is done in the simplest possible case for which gravity is described by 
general relativity and the matter content has a single component, namely a scalar field in a closed 
geometry. We show in particular that the spectrum of scalar perturbations can be affected by the 
bounce in a way that may depend on the wave number, even in the large scale limit. This may have 
important implications for string motivated models of the early Universe. 
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I. INTRODUCTION 

Although the idea that the universe could have expe- 
rienced a bounce in its remote past is old Q, it has 
recently come under new scrutiny 0, 0, lElSl with the 
advent of string motivated scenarios of the pre big bang 
kind 0, 13 ■ The main reason for this renewed interest 
is the fact that the most popular extensions of the stan- 
dard model of high energy physics, such as string or M 
theory, when applied to cosmology, i.e. in a four dimen- 
sional time dependent background, can lead to solutions 
with bouncing scale factors (see, e.g. Refs. 0, S and 
references therein). 

A crucial property to decide whether these new mod- 
els can be turned into realistic alternatives to the infla- 
tionary paradigm which, so far, has been so successful is 
the behavior of cosmological perturbations around these 
bouncing backgrounds. In particular, an important test 
is to calculate the evolution of the power spectrum of pri- 
mordial fluctuations through the bounce in order to see 
whether it can be made close to scale invariance, i.e. if 
there is any possibility, given the prebounce era, to get a 
Harrison-Zel'dovich power spectrum. 

From a technical point of view, the previous question 
is a nontrivial problem. Simple models, based on general 
relativity with flat spatial sections lead to the exis- 
tence of a curvature singularity at the bounce itself and 
therefore do not seem to represent viable physical models. 
In addition, it is difficult to understand how meaningful 
a perturbative scheme around a singular solution would 
be (see however Ref. 0), so we shall assume that the 
question of the calculation of the cosmological perturba- 
tions cannot be addressed in this way (see Ref. for 
more detailed discussions). In fact, it is believed that, in 
the vicinity of the bounce, string corrections become im- 
portant Typically, these corrections add to the grav- 
ity sector terms such as i?^, R^vpxR^"''^ , where RgpoX 
is the Riemann tensor and R the curvature scalar 
and higher order terms in the curvature. The effect of 
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these terms is, except in some specific instances 0, to 
smooth out the singularity ^3 (see also Ref. 0); this 
is, from a physical point of view, satisfactory and ex- 
pected. Unfortunately, one can show that the process 
of adding up more and more high curvature corrections 
does not lead to convergence towards a single solution, 
i.e. this solution explicitly depends on the choice of the 
stringy corrections [13 ■ This means that, at each order, 
the bouncing scale factor looks completely different from 
the scale factor obtained at the previous order. Never- 
theless, at a fixed order in the string corrections, one 
can in principle compute how the perturbations propa- 
gate through the bounce. The main disadvantage of the 
procedure is that it renders the computation extremely 
complicated and only numerical calculations are available 
to make the problem tractable. 

An important point to be noticed is that, as already 
mentioned above, most models assume that the spatial 
sections are flat all the time whereas, at the bounce, the 
curvature term is expected to play a crucial role. There- 
fore, it seems that for a bouncing universe, one cannot 
just throw away the curvature term because it does not 
play a significant role as it is the case for an inflation- 
ary universe. In fact, in that regard, the situation is the 
opposite of inflation: during the final stages of inflation, 
one can safely assume flat spatial sections because the 
three-curvature is getting more and more negligible as 
time passes, whereas even though the curvature may be 
negligible either in the remote past or in the future of 
the bounce, it has almost certainly no reason to be so in 
general. 

A way out to the previous difficulties, which would per- 
mit to undertake a tractable analytical calculation of the 
power spectrum, is the following. Far from the bounce, 
one usually considers the situation for which the curva- 
ture is small, even though the implementation of this 
particular point may not be in itself a trivial task. In 
this case, one can consider that, during the contracting 
and expanding phases, the spatial sections are essentially 
flat so that the well known results stemming from the 
theory of cosmological perturbations can be straightfor- 
wardly applied. Then, the main question becomes the 
effect of the bounce itself on the pre-bounce power spec- 
trum. Technically, this problem can be formulated as 
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follows 13 • Before the bounce, the perturbations are 
characterized by two modes, a dominant (denoted by D) 
and a sub-dominant {S) one. With k representing the co- 
moving wave number of a given Fourier mode, we write 
the k dependence of these modes as D^{k) and S-{k). 
In the same manner, after the bounce, one decomposes 
the perturbation as D+{k) and 5+(fc). The effect of the 
bounce is then entirely encoded into the form of the tran- 
sition matrix T(fc) defined by 



D+ 
S+ 



Til Ti2 
T21 T22 



(1) 



The mode of interest is of course the dominant mode 
in the expanding phase, D+{k) — Tii(fc)Z?_(fc) + 
Ti2{k)S-{k). This equation corresponds to a general "fc- 
mode mixing," i.e. the dominant mode after the bounce 
is a general linear combination of the dominant and sub- 
dominant modes before the bounce.^ 

A priori, various different situations can occur: the 
dominant mode in the contracting phase could acquire 
a scale invariant spectrum which is not conveyed to the 
dominant mode in the expanding phase because it turns 
out that Til = and T12 ^ ("fc-mode inversion," the 
scale invariant piece is passed to the "wrong mode" in 
the expanding phase); this is for instance what occurs if 
one applies the usual Israel junction conditions, known 
to apply for other cosmological transitions, at the bounce 
point 15]. 

Another possibility is that the dominant mode in the 
contracting phase be scale invariant but that this prop- 
erty is lost through the bounce due to a nontrivial k de- 
pendence of the coefficient Tn. Note that the opposite 
situation may also occur, for which the spectrum is ini- 
tially not scale invariant but is turned into it because of 
a nontrivial k dependence of the transition matrix. In 
fact, the common view concerning these last possibili- 
ties is that, for scales of astrophysical interest today, the 
bounce, lasting a short time, is expected to have no no- 
ticeable effect on those large scales Q . This is sometimes 
argued to come from general arguments such as "causal- 
ity," a point which is discussed thoroughly in Ref. [l^ . 
For instance, this is a basic assumption in the pertur- 
bation spectrum calculations in the pre big bang sce- 
nario • Technically, this means that the transfer matrix 
is assumed not to depend on k ^^l- Within this frame- 
work, the goal reduces to finding situations for which a 
scale invariant spectrum is produced in the contracting 
phase, and to ensure that this spectrum is passed to the 
dominant mode in the expanding phase, i.e. to insure 
that the matching conditions at the bounce do not imply 
a /c-mode inversion. 



^ We have introduced the notation "fc-mode mixing" in order to 
emphasize that this mode mixing is vahd for a fixed Fourier mode 
and should not be confused with the mode mixing coming from 
nonlinearities which involves modes with different wave numbers. 



The present article aims at examining whether the as- 
sumption that the transfer matrix is fc-independent is 
generically valid or not. For this purpose, we need to 
specify a class of models where bouncing solutions are 
possible and which allows simple analytical treatment of 
the perturbations through the bounce. We choose gen- 
eral relativity, positive curvature spatial section (see the 
remarks above), and describe the matter content by a 
scalar field; a similar strategy was used in Ref. |^. We 
do not assume anything relative to what happens away 
from the bounce, and in particular one could envisage 
that there the curvature is negligible; note that, as we 
show below, this implies the existence of a new phase. 
Therefore our closed geometry bounce can be viewed as 
an example of a transition connecting the contraction 
phase to the expanding phase with flat spatial sections, 
as already considered in the literature. 

This article is organized as follows. In the following 
section, we set the precise model and derive the basic 
equations both for the background and the perturba- 
tions. We then discuss how one can model a bounce 
in this framework and derive an explicit form for the po- 
tential of the scalar part of the classical perturbations 
(Bardeen potential) whose properties we then examine 
in details. This leads us to the main calculation of this 
article, namely that of the transfer matrix of Eq. We 
show that this matrix depends on /c in a nontrivial way 
provided that the null energy condition (NEC) is very 
close to being violated at the bounce. This illustrates, by 
means of a specific example, that the general argument 
according to which the limited but nonvanishing bounce 
duration could not affect the spectrum of long (i.e. longer 
that the duration itself) wavelength modes, is incorrect. 
We conclude by discussing this result, also showing that 
in the case under consideration, the propagation of gravi- 
tational waves (tensor modes) is qualitatively different of 
that of scalar modes since the former are never affected 
by the bounce. 



II. BASIC EQUATIONS 

We assume that the background model is given by 
a Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 
verse, i.e. 



ds^ =0^(77) 



1 -/Cr2 



-sm 



(2) 



In this equation, the constant parameter K, can always be 
rescaled such that JC — 0, ±1 and describes the curvature 
of the spatial sections. The time 77 is the conformal time 
related to the cosmic time by At — a{ri)dri. The matter 
is described by a homogeneous scalar field >p{ri) and the 
corresponding energy density and pressure respectively 
read 



/2 



2a2 



2a2 



V{ip). 



(3) 
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A prime denotes a derivative with respect to conformal 
time. The function V{(p) represents the potential of the 
scalar field. Einstein equations relate the scale factor to 
the energy density and pressure of the scalar field accord- 
ing to 



with the three-dimensional metric. The eigenfunctions 
fn{x^) of the Laplace-Beltrami operator on the spatial 
sections satisfy the equation 



-n{n + 2)/„ 



(8) 



:{2H'+H^ + IC) = KP, 



(4) 
(5) 



where we have defined 7i = a' /a and k = Sn/mpi, mpi 
being the Planck mass. The quantity p + p is then given 
by 

^,(p+p)^l^n^r^\{n^^n' + ic) = K^. (6) 

From the above equation, we see that the function r(77) is 
defined by T ee l-H' /H^+JC/H'^. It is directly related to 
the equation of state parameter lj = p/ phy the following 
relation w= (2r/3)(l -h /C/H^)"^ - 1. The function T (77) 
reduces to a constant for constant equation of state and 
is zero in the particular case of the de Sitter manifold. 
At the bounce, the Hubble parameter vanishes, TC = 0, 
while H' > 0, and therefore the only way to preserve the 
null energy condition p -|- p > is to have K. > 0. This 
is why, in this article, we restrict ourselves to the case 
K, > 0, i.e. the spatial sections are 3-spheres. Let us also 
notice that, being given a bouncing scale factor 0(77), it 
is sufficient to check that F > at all times in order for 
the scale factor to be a solution of the Einstein equations 
with a single real scalar field. 

A universe with closed spatial sections is characterized 
by two fundamental lengths. The first length is the Hub- 
ble length defined by £^ = a^/a' = a/a = (a dot 
denoting a derivative with respect to cosmic time t) and 
the second one is the curvature radius, = a/y/\JC\. 
The flat limit is recovered when £^ ^ as revealed 
by the equation |1 — f7| = £'^/£'^, where is the ratio 
of the total energy density p to the critical energy den- 
sity. When it comes to numerical applications, let us 
recall that one can safely assume the preferred value 0| 
Ho = lOO/i km-s~i-Mpc-i with h = O.lltom: leading 
to a Hubble distance scale now of ^ SOOOft.^^ Mpc ~ 
4.2 ± 0.2 Gpc. Moreover, with f^now = 1-02 ± 0.02, one 
has a curvature length, namely the scale factor as mea- 
sured now of order ao ^ 15h~^ Gpc (with K. ~ 1), 
the limit coming from the maximum allowed value for 



at one a level. 



At the perturbed level, and in the presence of density 
perturbations only, the metric takes the following form 



a2(77)|- (1 + 20) drf + 2d^BdTjdx' 



(1 - 24,) 7^^^ + 2V,d,E 



dx'dx^ 



}, (7) 



where is the metric of the spatial sections and the 
symbol denotes the covariant derivative associated 



where n is an integer. Note at this point that it is 
because of our normalization with a dimensionful scale 
factor a(r/), and hence dimensionless coordinates (77,0;*), 
implying a dimensionless operator itself, that the eigen- 
values of A are dimensionless integer numbers; with a 
different convention, i.e. with a dimensionless, one would 
have [A] = and an extra factor would appear in 
the right hand side of Eq. ijSJi. 

The modes n = 0, corresponding to a homogeneous 
deformation, and n — 1, being nothing but a global 
motion of the center of the 3-sphere, are pure gauge 
modes (Ts|: we will accordingly consider only values of 
n such that n > 1. In fact, for the relevant cosmological 
parameters discussed above, one finds that the values of 
n corresponding to characteristic distance scales of cos- 
mological interest now, namely lO^^h^^ Mpc < Dcosm S 
10'^h~^ Mpc, range between 30 and 3 x 10^ for the largest 
possible value of the total density now. For a reasonable 
value of f2now ~ 1.01, we find that n is between 60 and 
6 X 10^. 

The scalar perturbations are described by the four 
functions, (j), B, ip and E and, from them, it is possible 
to construct two gauge-invariant quantities, called the 
Bardeen potentials, and defined by (20li2lj| 



a{B - E') 



^--[B-E'). (9) 
a 



For simple form of matter with no anisotropic stress (this 
is the case for a scalar field), we have $ = 4'. Notice that 
the form of the Bardeen potentials is the same whatever 
the curvature of the spatial sections is. This is related to 
the fact that, even if /C > (or /C < 0), the FLRW metric 
remains conformally flat and the components of the per- 
turbed Weyl tensor remain unchanged. Then, Stewart 
lemma guarantees that the Bardeen potentials are still 
defined by the same equations (22l |. 

For the matter sector, the scalar field is written as tp-l- 
5tf{rj,x'^) where 6ip{ri,x'^) represents the inhomogeneous 
fluctuations. These fluctuations can be described by the 
gauge invariant quantity i5(^(si) = g^p ^ ^'{B — E'). 

The full set of Einstein equations can be written in 
terms of the gauge invariant quantities $ and only. 
Combining these equations permits to derive a master 
equation for the Bardeen potential (for (p' ^ 0) which 
reads 



<i>"+2{H 



<^'+ n{n+2)+2[n'-n^-2JCj $ . 

^ . 

This equation can be cast into a more convenient form. 
For this purpose, one introduces a new gauge-invariant 
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quantity, m, related to the Bardeen potential $ by 



where the function d is defined by 



(11) 



f 



a \p + p 



1/2 



f - 



Then, the equation of motion for the quantity u takes 
the form 



i(n + 2)---3/C(l-c2) 



In the above equation, one has 

"^ = 7 = "n^^ 



u = 0. 



(13) 



(14) 



for the scalar field, when use is made of the Klein-Gordon 
equation 



d(p 



= 0. 



(15) 



The quantity Cg of Eq. (|14|l can, in some regimes, be 
interpreted as the sound velocity. Let us now see how 
the flat case is recovered. The term &' jd is of order 7i^, 
namely d" jQ ~ This is a rigorous statement if 

the scale factor is a power law of the conformal time, 
which explains the usual confusion between the potential 
and the Hubble scale. However, this identification suffers 
from important exceptions, particularly relevant in the 
present context; see the discussion in Sec. lVCl Then the 
above equation can be re-written as 



47r" 1 
A2 ~ P~ ^2 (1 



_3_ 



u~ 0, 



(16) 



where we have used that the physical wavelength of a 
mode is A(ry) = 2'Ka{'q) / n{n + 2). In the limit ^ ^hj 
the last term of the equation becomes negligible. Then 
the equation of motion for u reduces to u" + [n{n -I- 2) — 
9" /9]u = where now 9 denotes the function defined 
previously in Eq. (|12|l but without the term proportional 
to K, (one can also show that, in the limit considered 
here, this term becomes negligible). Therefore, we have 
recovered the standard equation, valid for /C = 0. 

Note that in the flat limit for which /C = 0, the mode 
number n{n 2) that appears in Eq. H13I) appears to be 
a large number, and even more so when — > 1, as is 
the case in the usual cosmological calculations based on 
a period of inflation for which the approximation fc^ ^ 1 
is often done, permitting an expansion in powers of . 
This is not inconsistent though, because after N e-folds 
of inflation have happened, in a closed situation for in- 
stance, one expects i7 — 1 ^ c^^^, with > 55 to set 
the scales, and hence a gigantic value for the scale fac- 



tor normalization oq — Hq (il — 1) 



-1/2 



However, one 



then assumes, rightly, that the universe is almost flat, 
and chooses in general a different normalization for the 
scale factor, namely qq ~ Hq^, which is many orders of 
magnitude below what it ought to be. The relevant wave 
numbers, seen as eigenvalues of the Laplace-Beltrami op- 
erator, then must scale correspondingly, i.e. they are re- 
duced by the amount 



n{n + 2) 



(17) 



which is, indeed, much smaller than unity in any infla- 
tionary scenario. 

The equation of motion for the quantity u, Eq. H13|l . 
has the traditional form of a parametric oscillator equa- 
tion, i.e. of a "time-independent" Schrodinger equation. 
However, in the case /C 7^ 0, the effective potential cannot 
be written as the second time derivative of some function 
{9 in the flat case) over the same function because of the 
presence of the term — 3/C(l — c^) in the time-dependent 
frequency. This has for consequence that, on "super- 
horizon" scales where the term n{n + 2) is negligible, the 
solutions are not easily found contrary to the flat case 
where they are just u = 9 and u = 9 J <1t/9^ . 

So far, we have discussed the quantity u which is, up 
to some background functions, the Bardeen potential. In 
the framework of cosmological perturbations, there exists 
another important variable, usually denoted w, that we 
now consider. This quantity is important because its flat 
case equivalent naturally appears when one studies cos- 
mological perturbations of quantum-mechanical origin. 
In other words, this quantity is interesting for setting up 
physically well-motivated initial conditions whenever the 
curvature is negligible. Its definition reads 



3/C 



1 



5if 



(gi) 



H 



ICip' 



(18) 



n{n + 2) 



For /C = 0, it reduces to the well-known definition. The 
presence of the factor n{n + 2) in the definition above 
suggests however that this variable, in the /C ^ case, 
is not the canonical field that should be quantized to 
get initial conditions. This quantity involves (5(/3^s'-' and 
Since there are related by the perturbed Einstein 
equations, there is in fact only one degree of freedom as 
expected. The equation of motion for v reads 



n{n + 2) 



-3/C(l-c^) 



w = 0, 



where the quantity z is defined by 

_ a(p' 



Ha 1 - 3/C 



n{n + 2) 



(19) 



(20) 



This equation was obtained previously in Ref. [23 ■ The 
same remark as for the u equation applies: in the /C ^ 
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FIG. 1: Scale factors as functions of the conformal time rj 
corresponding to the de Sitter- hke solution [Eq. I)21|l . full 
line] and its various levels of approximations stemming from 
Eq. I30L namely up to quadratic (dashed), quartic (dotted), 
sixth (dot-dashed) and eighth power (dot-dot-dashed). The 
last two approximations, although clearly better from the 
point of view of the scale factor, do not lead to any new quali- 
tative information as far as the evolution of the perturbations 
is concerned. 



the effective potential is not only z" / z (as in the flat 
case) but is corrected by the +3IC{1 ~ c^) term. In ad- 
dition, the effective potential now depends on the wave 
number n through the quantity z. This has to be con- 
trasted with the effective potential for u, 0" /d, which is 
71-independent . 

Having defined the various quantities needed to study 
the evolution of cosmological perturbations through a 
bouncing phase, we now turn to the description of the 
bounce itself. 



III. MODELING THE BOUNCE 

In this section, we define precisely the behavior of the 
scale factor during the bouncing epoch, then discuss its 
relation with the following eras of standard cosmology 
and derive the relevant perturbation potentials. 



A. The de Sitter-like bounce 

Once the background is fixed, the effective potentials 
for the quantities u and v are completely specified. In 
this section, our aim is therefore to discuss how one can 
model the scale factor of a bouncing universe. At this 
point, one should notice the differences (and similarities) 
with inflation. In an inflationary universe, the behavior 
of the scale factor is known: essentially, this is a oc |ry|~^, 
i.e. the de Sitter phase. However, one can also treat 



FIG. 2: Behavior of the scalar field and its coordinate time 
derivative as functions of the conformal time rj (varying be- 
tween —n/2 and n/2 for the overall evolution of the Universe) 
for the de Sitter-like solution with 770 = 1.01. 



slightly more complicated backgrounds by means of an 
expansion around this de Sitter solution. This expansion 
is characterized by the so-called slow-roll parameters , 
which are constrained to be small. The de Sitter solution 
also exists in the bounce case ^3 ^^.d, as we shall see, it 
can be used in much the same way. However, contrary to 
the inflation case, there is no fundamental reason why the 
background equation of state should be close to vacuum. 
Despite this fact, one can nevertheless expand around 
the /C = 1 de Sitter spacetime and similarly deflne pa- 
rameters which control the departure from it. Obviously, 
those parameters are not subject to tight constraints, and 
in particular are not required to be small. 

For /C > 0, the de Sitter solution corresponds to 
the scale factor a{t) = aq cosh(u;t), which is expressed as 
a function of the cosmic time i, with lu = l/ag. More gen- 
eral solutions are obtained by relaxing this last constraint 
and considering a general value for u. These de Sitter- 
like solutions are the ones we shall be concerned with in 
what follows: our expansion will be based on these solu- 
tions. In terms of conformal time, one can integrate the 
relation adr] = dt to get 



0(77) = aoyi + tan2(^-^^, (21) 

where the conformal time is bounded within the range 
— 7r/2 < 77/770 < 7r/2 and the conformal time duration 
770 is related to the de Sitter coefficient uj through 770 = 
{aoUj)~^ [the solution is shown in Fig.p. 

In order to understand the dynamics of this solution, 
one needs to obtain the evolution of the scalar field. It 
can be integrated straightforwardly with the scale factor 
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FIG. 3: The shape 1271 of the potential for the scalar field ip 
(in units of the Planck mass = mp|/\/87r) for different 

values of the bounce characteristic conformal time 770 • The full 
lines are respectively for t^o = 1.001 (above) and rjo — 1.01 
(below), the dashed line corresponds to r;o ~ 1-1, and the 
dotted line is for 770 = 1-5. In the strict de Sitter limit -qo —> 1, 
the potential goes to the constant value V{(fi) = 3/(Kao), 
which explains why the 770 = 1.001 seems almost constant as 
it oscillates with a very small amplitude around its central 
value [(3 — l/?7o) in these units]. 



(l?T|l : from Eqs. Q and one obtains 



/2T 

^ = fa + y — [V + ^"^0 



(22) 



where we have set (p —^ (po as the cosmic time t — > —00, 
i.e. as 77/770 — Ti'/S, and we also have defined a param- 
eter 



Tee1-4 



(23) 



for further convenience. We shaU keep this definition 
later on for more general bounces than the quasi-de Sitter 
ones. 

It should be noted that the parameter T, in the case of 
de Sitter like expansion (|21|l is, according to the definition 
©, Tds = H^r, which is proportional to p + p. As a 
result, the null energy condition at the bounce can only 
be satisfied provided T > 0, i.e. if [770 1 > 1: indeed, one 
has 



lini (p 



T 

P) = 2^, 
"0 



(24) 



a relation which we shall use in the rest of the paper to 
define T in a solution-independent way. As emphasized 
before, the case 770 — 1 corresponds to a constant scalar 
field potential and to an equation of state p — ~p and is 
thus the exact counterpart of the inflationary de Sitter 



solution. The scalar field time derivative is now simply 
obtained as 



1/2 



lit 



dip 
adi] 



1 



2T 



1 + tan^ 



(25) 



Both the field and its time derivative are displayed in 
Fig. [3 for a case having 770 ^ 1 as functions of the con- 
formal time. 

It is now a simple matter to derive the corresponding 
potential for the scalar field which solves Einstein equa- 
tions Q and It reads 



KalV{ip) 



n' + 2 {n^ + jc) 



(a/ao) 

i.e. , with the solution (|?T|l above. 



KalV{(p) 



2T sin^ 



v/2k 



T-V2(^_^o) 



(26) 



, (27) 



and it is displayed in Fig. |31 with a specific choice of ini- 
tial conditions for the field. From Figs. |21 and O one 
sees that the universe starts at either a maximum or a 
minimum of the potential, in both cases with a non- 
vanishing amount of kinetic energy in the scalar field 

(^'(-^770/2) - V2t7^- 

One remarkable property of the above model is that 
the effective potential for density perturbations remains 
very simple even if 770 7^ 1. Indeed, assuming the de Sitter 
like solution (|21l) and plugging it into the form (|12|l yields 



ao 



2(%^-l) 



(28) 



which, together with — —1/3 [this stems from Eq. H14() 
with the solution (|21|l and the scalar field (|22|) ] leads to 



9 ' 



(29) 



or, in other words, the potential for the variable u does 
not depend on time for the de Sitter-like solution. Be- 
sides, the maximum value achievable by this potential is 
given for the limiting case 770 — > 1, and it is y^'^'^'™-^'^^ — 3; 
this potential can only interact with the modes 77, = and 
71 = 1, which we already mentioned are gauge modes. 
This was to be expected since the exact de Sitter solu- 
tion, in this bouncing situation as in the more usual in- 
fiationary scenario, does not amplify scalar perturbations 
by any amount. We believe that this is what happens in 
Fig. 4 of Ref. Q in which the Bardeen potential mode 
71 = 10 is seen to oscillate while passing through a de 
Sitter-like bounce, refiecting the nondomination of the 
potential at this point. 
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B. General bouncing scale factor 

We now assume that the universe experiences a regu- 
lar bounce at the time 77 = 0. This means there exists 
a particular function a(r/) which can always be Taylor 
expanded in the vicinity of r/ = 0. Since we are inter- 
ested in understanding the behavior of the perturbations 
through the bounce, and because the effective potentials 
for density perturbations involve derivatives of the scale 
factor only up to the fourth order, a description of the 
scale factor up to rj"^ only is sufficient. We therefore set 



1 



(30) 

which defines the parameters ao, the radius of the uni- 
verse at the bounce, rjo, the typical conformal time scale 
of the bounce, 6 and ^. They control the magnitude of 
each term of the expansion. 

For the scale factor ()30|l to be a solution of Einstein 
equations with a scalar field as matter content, the func- 
tion 0(77) must be chosen such that tjq is greater than 
unity, which is only a necessary condition as discussed 
below. The parameters oq and rjo also provide the tan- 
gent de Sitter-like solution H21|) . whereas 6 and ^ mea- 
sure the deviation with respect to this de Sitter-like so- 
lution. Eq. H30|l represents a double expansion, both 
around 77 = and around the de Sitter-like solution dis- 
cussed in the previous section since 5 = ^ = exactly 
corresponds to the small 77 expansion of the scale factor 
(|21|l : this explains, among others, the factor 5/24 in this 
equation. The parameters 6 and ^ are in a certain sense 
similar to the traditional slow-roll parameters. Both the 
de Sitter form and its approximations are shown in Fig.^ 
as functions of the conformal time. 

We now discuss how the general expansion H3U|) can 
be related to an underlying particle physics model. The 
scale factor is entirely specified once the scalar field po- 
tential V{lp) and the initial conditions, for instance the 
values of the scalar field and its first derivative at the 
bounce: </?o, fo, have been chosen. This means that 
there exists a relation between these last quantities and 
the parameters oq, ?7o, S and ^ characterizing the expan- 
sion H3(J|I. We now establish what this relation is. This 
can be done easily by solving the Einstein equations Q 
and lO in the vicinity of the bounce. In practice, we 
insert the Taylor expansion of the scalar field, 

ip{r])^ipo + ip'oV+l^W + Oiv^), (31) 



and of the scale factor, Eq. in Eqs. Q and © and 
identify, order by order, the various terms appearing in 
the resulting expressions. To zeroth order, this gives 



6 



1 



The last relation can also be re-written as 
T 



12 



(32) 



(33) 



as expected from Eqs. ^ and H24|l . We see that the mag- 
nitude of the scalar field conformal time gradient at the 
bounce determines the value of the parameter T. Typ- 
ically, one expects T ^ 1 since the order of magnitude 
of the scalar field and its derivatives should be such that 
V?o <C TTZpi in order for the field theory to make sense. 
If the velocity of the field vanishes at the bounce, then 
T = 0. 

To first order, the Einstein equations yield 



AV 



2dy 



12(5 



0, 
0. 



(34) 
(35) 



In the following, we will be mainly interested in the sit- 
uation where the bounce is symmetric, that is, we shall 
demand that (5 = 0. Then, there are two ways of satis- 
fying the Einstein equations. Either the kinetic energy 
vanishes at the bounce or i/Jq 7^ but then (^g ~ and 



0. This means that the bounce occurs at the 



minimum of the scalar field potential. This also implies 
that, in this case, the minimum of the potential cannot 
vanish, T^(</3o) 7^ 0, see Eq. If, for instance, the po- 

tential is given by y{}f>) oc (^", as is the case for instance 
of the model studied in Ref. , then the only way to have 
a symmetric bounce is to satisfy the condition (pg = at 
the bounce and, as a consequence, one necessarily has 
T = 0. In the following, we will be mainly interested in 
the second situation, i.e. ipQ ^ 0, since we will show that 
amplitude of the spectrum is controlled by the parame- 
ter T. In this case, going to the next order allows us to 
determine what the parameter ^ is. The result reads 



1 



5(2-K(^(2)2,cF2(^g) 



(2 - Kip'^)Vo 



,'2\2 



(6 - K(p'o ) 



-(6 - >iip'o' f^o 



(36) 



Let us notice that this parameter depends on the second 
order derivative of the potential at the bounce. 

Assuming a symmetric bounce for now on, i.e. setting 
6 = 0, some restrictions can be put on the numerical 
value of ^. They stem from the fact that we demand 0(7;) 
to be positive in the range —7^0 < rj < rjo and to describe 
a bounce, i.e. a' > for < 7; < 770. This latter condition 
turns out to be more stringent and implies ^ > —11/5. 
Moreover, if we further require that the scale factor H30() 
be solution of Einstein equation sourced by a single scalar 
field, we see from Eq. © that H^T must be positive. 
Around the bounce, this is 



^-2^(1 



2T) ( ^ 



O 



(37) 



which will be positive definite in a small but finite neigh- 
borhood of 77 = provided ^ < in the limit T —s- we 
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will be concerned with. Combining both constraints, we 
arrive at 

-y < e < 0. (38) 

The approximation method we discuss later does not al- 
low to consider very small values for ^, so that in practice, 
we shall use -11/5 < ^ < -0.1. 

We now assume the fiducial expansion H3U|) for the 
bounce through which we want to propagate the per- 
turbations. Let us however first examine the connection 
of this bounce to the standard cosmological epochs of 
radiation and matter domination. 



between the bounce and the standard hot big bang, an- 
other phase must necessary occur whose main effect will 
be to drive TC to sufhciently large values. This is usually 
the role played by a phase of inflation. 

With the general framework thus clarified, let us turn 
to the evolution of the scalar gravitational perturbations 
through the bounce by means of evaluating the effective 
potential for the variable u related with the Bardeen po- 
tential through Eq. (|ll|l . We discuss the potential for the 
variable v in the discussion section IV Al below. 

D. The potential Vu{ri) 



C. Connecting the bounce to standard cosmology 

In this section, we study how the bounce that we de- 
scribed previously can be connected to an epoch of the 
standard hot big bang model. In particular, we study the 
connection with a radiation dominated era. In this case, 
the scale factor can be written as 



a{ri) = flr sin(?7 — 77r) 



(39) 



where and yyi- are two parameters to be fixed with the 
help of the matching conditions. We match this scale 
factor to the bouncing scale factor given by Eq. H30|l . us- 
ing the junction conditions, known to be valid even in 
the curved spatial section case, as derived in Ref. [2fil |. 
namely [a] = [a'] = 0. The matching is performed at 
r] = Tjj such that ryj <^ r/o in order for our quartic approx- 
imation of the scale factor to still be meaningful. The 
matching conditions imply that the Hubble parameter at 
the matching time is given by 



1 + -x^ 



(40) 



where x = <C 1. From the above formula, one 

sees that it is not possible to connect the bounce to an 
epoch where the curvature is negligible, provided the null 
energy condition, which demands rj^ > 1 [see discussion 
around Eq. H24(l ]. is still satisfied. As a consequence, 
this implies that Ti-irj-j) cannot be large in comparison 
to /C 1; in fact, since t^q — 1 and a; 1, Ti.^ is ex- 
pected to be negligibly small compared to unity right 
after the bounce. This means that one necessarily con- 
nects the bounce to a regime where the curvature is im- 
portant or, in other words, in a region where the sine 
function appearing in the scale factor H39|l cannot be ap- 
proximated by the first term of the Taylor expansion, 
a (77) ~ a-[{rj — rj-n) . The only way to avoid this conclusion 
would be to violate the null energy condition, as already 
noticed in Ref. and to have a small 770 but then it 
would have been useless to consider the case /C = 1 for 
modeling the bounce since this was done precisely in or- 
der to satisfy this condition. Therefore, we conclude that 



The effective potential for the variable u in the de 
Sitter- like solution is, according to Eq. 1)29(1 . constant in 
time. This is however very specific to this particular so- 
lution, as any displacement away from it immediately 
leads to a different form of the potential. This is illus- 
trated in Fig. 21 which shows the relative accuracy of the 
expansion H30() around the de Sitter- like solution H21|l . 
It is also clear from the figure that the expansion H30|l . 
if pushed to sufhciently high orders in 77, gives back the 
correct constant value over a large range of conformal 
times. Let us now turn to the more general bounce case 
of Eq. jUJ. 

Arbitrary values for the parameter ^ restricted to the 
range of interest discussed above lead to the generic shape 
illustrated in Fig. [S| The calculation of the effective po- 
tential is extremely complicated even with the quartic 
approximation of the scale factor. Even if it can be done 
in full generality since, for a scale factor given by Eq. I|30|l . 
the potential Vu{ri) reads 



3/C (1 - c^) = 



Q24{V) 



(41) 



where P2i{'n) and Q2A{'>i) are two polynomials of order 
24, in practice the calculation is not tractable. However, 
since in practice we always have 77/770 1; only the first 
monomials are important. One can check that the fol- 
lowing approximation 



Co + C2'q 



do + ^2??^ + diV 



(42) 



is extremely good, see Fig.|Hl In this expression, we have 
only kept the first two monomials at the numerator and 
the first three ones at the denominator. The coefficients 
c, 's and d, 's can be written as 



Co 

C2 

do 
d2 
di 



= 2772(772 - 1)(2 - 10772 + 8774 + 50 



(43) 



49 + 48776 + 55$ + 50^2 - 1877^(6 + 50 
+772(11 + 350, (44) 
12770^(772 - 1)2 (45) 
12772(772 -l)(-3 + 3772 -50 (46) 
-1 + 1577^^(4 + + 5e(16 + 150 - ^/o(59 + 950 • 

(47) 
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FIG. 4: Absolute value of the effective potential Vu{ii) for the perturbation variable u{ri) for the de Sitter- like case (full line 
on both panels), for which it is constant and for the various approximation levels (from quadratic to eighth power of the scale 
factor). The left panel shows the potential as obtained by using the quadratic (dotted line) and quartic (dashed) expansions 
of the scale factor only, whereas the right panel presents the situation when quartic (dashed), sixth (dotted) and eighth (dot- 
dashed) terms are used. It is clear that the quadratic approximation is qualitatively wrong and cannot be used to describe a 
de Sitter bounce. The value 770 = 1.01 has been used to derive these plots. 
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FIG. 5: Absolute value of the potential Vuijf) as a function 
of rescaled conformal time rj/rjo for jyo = 1.01 as derived us- 
ing either the assumption that the scale factor behaves as a 

square root, i.e. a — ao^Jl + (J]/rjo)^ , (full line) or Eq. I30II 
up to quadratic (dotted line) and quartic order with 5 — 
and ^ = —2/5 (dashed line). The quartic approximation is 
extremely close to the exact solution, exemplifying its accu- 
racy, while the quadratic approximation appears to be at best 
qualitatively correct. 
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FIG. 6: The potential K(?;) [fuU line, Eq. (|4TJ] and its ap- 
proximation stemming from Eq. H42^ Vu'^'^'^\ri) (dot-dashed 
line) in full details with the same parameters as in Fig |3 
Also shown is the parabolic approximation that will be used 
in Sec. II V Al (dashed line), the various conformal times in- 
volved in the calculations in the text, as well as a pictorial 
definition of the regions I, II and III used for the matching 
of the perturbations also in Sec. IIV Al The mode n{n + 2) 
interacts, in this example, only with the central part of the 
potential. 



Vu{r]) ai r] = is given by 



Equipped with Eq. (|42|l . we can now compute the height 
and the position of the central peak and of the wings. 
Let us start with the central peak. The absolute value of 



3co ^ 2 - lOr;^ + 8tjI + 5^ 



(48) 



The most important property of the above formula is 
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that it diverges as ryo ^ 1- It is shown in Ref. that 
this property also holds in the case 5^0, and is there- 
fore generic, i.e. not restricted to symmetric bounces. We 
have seen previously that the values of n of astrophysi- 
cal interest are such that n ^ 1. Therefore, a necessary 
condition for the bounce to affect the spectrum of the 
fluctuations is that ryo be close to one. As already dis- 
cussed, the physical interpretation is that one must be 
very close to a violation of the null energy condition. In 
this case, it is more convenient to work with the variable 
T introduced in Eq. In practice, T must be a tiny 

number in order to get a modification of the spectrum. 
We have seen that this is to be expected since T is the 
square of the ratio of the scalar field conformal time gra- 
dient at the bounce to the Planck mass. The amplitude 
of T controls the maximum value of n below which the 
perturbation modes will be affected by the bounce. The 
crucial point is that for T <C 1, large scales, having cos- 
mological and astrophysical relevance, can be modified 
as they evolve through the bounce. 

Assuming for now on that T <C 1, it is sufficient to 
Taylor expand everything in terms of this parameter to 
get an accurate approximation. For Vq, one gets 



^0 



2T 



(3-50 + O(T). 



(49) 



Another interesting quantity is the time rj.^ for which 
Vu ~ 0, see Fig. El This time is given by -q^ = —ct^^jci 
which leads to 




(50) 



The above equations means that, in the limit T — *■ 0, 
the width of the potential goes to zero while its height 
increases unboundedly. Finally, let us describe the wings 
of the potential. The position of the wings can be derived 
from the condition = (77 7^ 0). This gives 



1 



C2d4 



and the Taylor expansion in T reads 



-^ + 0(T3/2). 
5^ 



(51) 



(52) 



One sees that 77^ — 2775. at first order in T. Therefore, 
77w also goes to zero when T tends to zero. The height of 
the wings is just given by ^^(ryw) and can be expressed 
as 



6T 



3+-0 +0(T). 



(53) 



The height of the wing also diverges as one approaches 
the violation of the null energy condition and, at first 
order in T, one has Vq/Vw — 3. This concludes the de- 
scription of the perturbation potential with which we now 
examine the fate of the perturbations themselves. 



IV. CALCULATION OF THE TRANSFER 
MATRIX 

The purpose of this section, which is also our main re- 
sult, is to show that the transfer matrix T of Eq. ^ may 
depend on the wave number n in a way which we derive. 
We found that two completely different and independent 
methods, one based on a piecewise expansion of the po- 
tential and the other assuming the potential to behave 
mathematically as a distribution rather than a simple 
function in the limit T — > 0, lead to comparable results, 
both for the final spectrum itself and its magnitude. We 
examine these methods in turn. 



A. Method I: piecewise solution approach 

When trying to evaluate the transfer matrix of Eq. 
with the potential Vu{r]) derived in the previous section, 
one immediately faces a difficulty, namely that unfortu- 
nately, even with the simple form given by Eq. (|42l) . the 
equation of motion of the variable u is not integrable an- 
alytically. However, one can find piecewise solutions. For 
a mode which interacts only with the central peak of the 
barrier, i.e. which is above the wings, the potential is es- 
sentially zero for \ri\ > rj^,. This corresponds to regions I 
and III in Fig.|Hl In the central region, region II in Fig.O 
corresponding to \r]\ < rj^,, we model the bounce by a 
parabola with a minimum at — Vb and which vanishes at 
77 — ±ri2. To summarize, our piecewise potential is given 

by 



-Vz<V<Vz, (54) 



In each region, the function u is the sum of two modes 
and can be expressed as 

Ui{n,r]) = Ai{n)fi{n,ri) + B^{n)g^{n,r]) , i = I, II, III. 

(55) 

Before and after the interaction with the barrier, the so- 
lution are plane waves. 




1 - 


ii)) 







1 



1 



Jkri 



2k 



(56) 



where we have introduced the quantity k = ^n(n 2) 
and the normalization is chosen such as to simplify fur- 
ther calculations (unit Wronskian). In region II, one has 
an even and an odd mode, i.e. 



/ii {-n) = fu iv) , 9u i^v) = -9u iv) ■ 



(57) 



For the moment, we do not specify what fuirj) and guirj) 
are since we are trying to keep the calculation as gen- 
eral as possible (for example, we could imagine other 
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parametrization of the potential in the central region for u and its derivative u' a,t rj = ztrj^. Straightforward cal- 
which /jj and would be different). Our goal is to pre- culations lead to 
diet what ^jjj and _Bjjj are. For this purpose, we match 



An 



2ikri^ 



2ikW{n) 



I 





'A' 







, (58) 



where W{n) is the Wronskian of the function /jj and 
namely W{n) — fng[-^ — f^-^gu- In the previous expres- 
sions, all the functions are expressed at the point r] = rjz 
(we have used the parity properties of the function /,j and 
(jjj). The above matrix is general and is parametrized by 
only four numbers: /„, /^^ and g[^. Any model per- 
mitting to calculate what these numbers are allows us to 
estimate the transfer matrix on the bounce given above. 

We now use the parabohc model introduced before. 
If we perform the following change of variable, 77 = 
\/ry77(2^7^yx, then the equation of motion for u in re- 
gion II takes the form 



d^~ T + "'"^'^' 



where the parameter a is given by 



1 



1 



n{n + 2) 



(59) 



(60) 



Equation (|59|) can be solved exactly in terms of cylinder 
parabolic functions |27j . Since the potential is symmet- 



ric, the solutions can always be chosen to be even and 
odd. The explicit expression of the even and odd solu- 
tions are respectively 



/„(^) - e-^/^''^/(^-),F,(| + i;i;^,^),(61) 
5„(,) = ,e-v^''^/(^-)iF,(| + J;|^,^)(62) 



where i^Fi is the Kummer confluent hypergeometric func- 
tion. As already mentioned previously, these functions 
and their derivatives must be evaluated at rj = rj^ and 
then expanded in the parameter T. The first step is to 
calculate the parameter a. This gives 



^ ^ L + {-61 + 5[-53 + 8n(n + 2)]a ^ ^ ^^^^^ 



2V2 



200V2^2 



(63) 

Using this expansion and that of Vq and rjz, one obtains 
at first order in T 



/ll (??z 



\.^]+ 0{T'/') ^ 0.798 + 0(Ti/2) , 



V- 


-5^ 


V- 


-5^ 



-V21F1 



_ojiiV5| 

Ti/2 

1 -1/(2V2) 



0(T^/^), 
(6-3V2)iFi 



2-V2 1 1 
8 ' 2' V2 



6-\/2 3 1 



2-1)1^^1 



(3V2~l)iFi 



10-^/2 3 1 
8 '2' 71 



14 - \/2 5 1 
8 '2'71 



Ti/2 



(64) 
(65) 

0(Ti/2), (66) 
(67) 



0.878 + 0(Ti/2). 



(69) 



r 



The expression for the derivatives can be easily recovered of a Kummer hypergeometric function, iF[{a,(3,z) = 
if one uses the following expression giving the derivative (a//3)iFi (a -|- 1,13 + 1, -z), where a prime in this context 
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means a derivative with respect to the argument z of the 
hypergeometric function. 

The next step consists in inserting these relations into 
the general form of the transfer matrix and then in ex- 
panding the resulting expression in the parameter T. The 
result reads 



-0.624i. 



1 1 



n(n + 2) -1/ Ti/2 



1 



(70) 



Several remarks are in order at this point. First, the 
formula above applies only for the modes actually inter- 
acting with the potential, namely those having n(n -\- 
2) < Vb, otherwise, T-u is obviously the identity. Note 
also that, in the former case, the transfer matrix is 
n-dependent. This means that the bounce affects the 
spectrum and therefore disproves a priori any general 
argument stating that the spectrum should propagate 
through the bounce without being modified 0, uM ■ 
sides, we see that the amplitude is divergent as T goes 
to zero. However, one should remember that we are not 
interested in the spectrum of u itself but rather in the 
spectrum of the Bardeen potential <&. The relation be- 
tween $ and u, Eq. (fTTI) . together with Eq. if^ . leads to 
the remarkable result that the terms in T cancel out ex- 
actly and that the resulting spectrum is T-independent, 
and thus perfectly finite even in the T — > limit. Fi- 
nally, the ^ dependence of the overall amplitude is also 
predicted by this calculation. As expected, there is no 
net effect in the limit ^ — at which the bounce is ef- 
fectively de Sitter and thus can amplify no amount of 
perturbation. In this last case, the calculation leading to 
Eq. H70|) is not accurate enough and should be done at 
a higher order in T since the leading order vanishes; one 
should then find that the transition matrix is essentially 
the identity (de Sitter) plus some correction vanishing in 
the hmit T ^ 0. 



B. Method II: Distributional approach 

We show in this section that the previous result can 
be understood in very simple terms and that the result 
of the previous section can be reproduced by a back-of- 
the-envelope calculation. The crucial observation is that 
the height of the potential Vu diverges as T goes to zero 
while its width shrinks to zero. This suggests that there 
is something like a Dirac i5-function at play.^ To study 



^ Note that even though we send T — > in this section, this is 
merely a computational artifact allowing an easy calculation of 
the effect. The true value of T must be nonvanishing, although 
tiny, so the calculation of this section is accurate only for those 
modes interacting with the potential. Therefore, the presence 
of the Dirac distribution in no way implies the existence of a 
singular behavior either of the potential Vu{ri) or of the modes 
u(rj) themselves as long as T 7^ 0. 



this point we calculate the integral of the potential. One 

gets 

/+00 r+vo 
[K(r)-4]dr ^ / V;(-PP)(r)clT (71) 
-00 J —rjv, 

+ C'(T°).(72) 



(73) 



8T 

Thus, the potential can be re-written as 
K(r?) = -CTAT(r?), 



where the constant Cy is given by Cx = [— Stt^^/ (BT)]-'^/^ 
and where the function Ax('7) is a representation of the 
Dirac (5-function, i.e. 



lim Axff?) = 5{ri) . 



(74) 



In a certain sense, the potential Vu{ri) possesses diver- 
gences "worst" than a Dirac (S-function. The equation of 
motion of the quantity u can now be written as 



u" + [n{n -f 2) -I- Cx(5(?/)]u = , 



(75) 



i.e. a well-known equation in the context of quantum 
mechanics. The matching conditions are [u] = and 
[u'] = — Cxu(O), the last one coming from an integra- 
tion of the equation of motion across a thin shell around 
77 = 0. This reduces to 



A,,, B,,, = A, — B, — 



Cx 



(76) 

:(4+i?J.(77) 



iy''n{n + 2 

Straightforward algebraic manipulations lead to the fol- 
lowing transfer matrix, under the assumption that the 
second term of the last equation dominates over the first 
since Cx ^ cxd as T — > 0, 



32n{n + 2) 



1 



1 1 



1 -1 / Ti/2 



(78) 



It is interesting to compare Eq. H70|) with Eq. (TSJ. 
The numerical coefficient in the above equation is 
7r\/5/(4V2) — 1.242, to be compared with the coefficient 
0.624 found in Eq. (|70|l . The difference is approxima- 
tively a factor 1/2 in the amplitude. This difference can 
be interpreted in the following way. When the matrix 
transfer is computed using the matching procedure, one 
uses the parabola formula for the potential and one ne- 
glects the wings of the potential. The area of the central 
part of the potential is given by 



K(T)dr = 



looe 4V2 



45T 



yi^pp)(T)d7 



(79) 

We see that there is factor 4\/2/ (Stt) between the area be- 
low the central part and the area below the whole poten- 
tial including the wings. Since the matching procedure is 
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sensitive to the central part only whereas the calculation 
of the Dirac iJ-function is sensitive to the whole potential, 
we therefore expect a factor 4-\/2/ (37r) between the cor- 
responding two amplitudes. We have 4-\/2/(37r) ~ 0.600 
and hence we recover approximatively the factor 1/2 
mentioned above. The correct amplitude is the one given 
by the Dirac 5- function calculation and is ~ 1.25. 



V. DISCUSSION 

We now complete the description of the propagation 
of perturbations through a general relativistic bounce by 
some considerations regarding the variables u and v, the 
spectrum of tensor modes and a comparison with other 
known transitions in cosmology. 



A. u versus v 

An interesting issue, debated at length in the litera- 
ture 0, 0, IE 01 7 is how the variables u and v behave as 
they go through the bounce. As a first step towards un- 
derstanding what is the variable that is the most useful, 
let us construct the potentials for both, as in Fig. [7| It 
is clear from this figure that the terms appearing in the 
potential for the variable v, namely 



= — -f 3/C(l-c2 



(80) 



never compensate each others, as was found to be the 
case for Vu of Eq. (|^ . so the resulting potential is di- 
vergent at some points which, furthermore, depend on 
the wavelength index n. This provides a first hint that 
V is not the correct variable to work with, and indicates 
that as one approaches the bounce, or as the curvature 
becomes non- negligible, v ceases to be the good quantum 
variable (see e.g. Ref. |2ll|). 

In the present context, it is easy to show that u and v 
are related by the following relation: 



1 



3/C 



1 



aVT ' 



(81) 



n{n + 2) 



We know from the previous considerations that u is con- 
tinuous and that u' may have a finite jump at 77 = pro- 
vided T is small enough. From the above equation, we 
conclude that the variable v possesses divergences during 
the bounce. These divergences are given by the zeros of 
the argument of the square root at the denominator of 
the previous equation. In other words, v diverges when 



= 1 



n{n + 2) 



(82) 




1.0 



FIG. 7: Construction of he potentials Vu and 14 for the per- 
turbation variables u and v in the special case of the square 
root form for the scale factor as in Fig. El [or Eq. ^ with 
5 = and ^ = —2/5, the corresponding curves being visually 
undistinguishable] with rjQ = 1.01. According to Eqs. 
and 119II . the potentials depend on three possible terms, 
namely 9" /O, SIC (l — Cg) and z" /z, respectively plotted as 
the dashed, dotted, and dot-dashed curves. The potential for 
u is also shown as the full line. The pole at 77 = in either 
6" /6 and SJC (l — c^) appears with opposite sign but is oth- 
erwise the same, so there is an exact compensation, so that 
the full potential is everywhere well-behaved. This is clearly 
not possible for the potential Vv, since there are more poles in 
z" /z than there are in SIC (l — Cg), so no compensation can 
occur at these points, but also the pole at 77 = appear with 
the same sign; the potential for the function v follows z" /z, 
up to small corrections and was therefore not plotted here. 



as found in Fig. [7| 



Some remarks are in order at that point. First, an in- 
teresting feature is that w = at the bounce (hence is 
regular) and that the divergences occur before and af- 
ter the bounce but not at the bounce itself even though 
its potential actually diverges at this point. This is 
because the effective velocity of sound diverges at the 
bounce. Secondly, the time at which v divergences is n- 
dependent as can clearly be seen from Eq. . Thirdly, 
the physical interpretation of this divergence is subtle. If 
V had the usual interpretation (i.e. the variable that is 
canonically quantized), the divergence would clearly be 
a problem. Roughly speaking, this would mean explo- 
sive particle creations and, as a consequence that there 
is a back-reaction problem. More seriously, this diver- 
gence would be at odd with the fact that the Bardeen 
potential remains finite and small. As already discussed 
below Eq. I|18|) . there are reasons to believe that, in the 
case K. — 1, the variable v introduced before is not the 
variable that appears in the action for cosmological per- 
turbations. This last variable should remain finite during 
the bounce. 
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B. Density perturbations versus gravitational 

waves 



words, we study this transition with /C = 0. The scale 
factor can be given the form 



The evolution of the tensorial modes of perturbations 
fjL = ah, where h is, roughly speaking, the amplitude of 
the gravitational wave, stems from the relation [2^l2^ 



A* 



n(n + 2)-/C 



M = 0, 



(83) 



i.e. an equation similar to that valid for the scalar modes 
but with a potential simply given by Vh = /C + a" /a. 
Within the framework of our bouncing solution, this is 



2 V'70 



2 V%/ 24 \-qo 



1. (84) 



which can be simply analyzed as follows. 

As direct calculation reveals, the potential Vh of 
Eq. (|84|l has either a single maximum located at = 
if the expansion parameter ^ < —4/5, or a minimum at 
77 = and two maxima at the points ?7max given by 



2 ^2 , ^5(5 + 60-1 



(85) 



provided —6/5 < £, < —4/5. In the latter case, the max- 
imum value attained by the gravitational wave potential 
is 



15(1 + 



2 + ^5 (5 + 60 



1, 



5 - ^ - 5' 



otherwise. 



Since 1^ < and 770 > 1, this means that the maximum 
value for the potential is less than 22/7 ~ 3.14 in all cases 
of physical interest. In other words, and since in these 
units the cosmologically relevant modes are those having 
n ^ 1, the potential is dominated at all times during the 
bounce itself, and therefore cannot lead to tensor mode 
production. There is therefore a qualitative difference 
between the tensor and the scalar modes since the latter 
can be affected by the bounce provided the NEC is almost 
violated, while the former are never affected, regardless 
of the underlying parameter values. 



C. Comparison with other transitions 

In order to make a comparison of our bouncing era 
to other known transitions, we first consider below the 
radiation to matter transition under the hypothesis that 
this occurs, as observation demands, at some time 77oq 
such that the three-space curvature is negligible. In other 



a{r]) 



\VcqJ 



(87) 



where 61 = 62 = V2— 1 is chosen such that a{r] = rjcq) — 
Acq- We have emphasized the two different normaliza- 
tion factors bi and 62 because they play a different role 
in the potentials for either u or v. Indeed, the poten- 

tial for . in this case is (aVr)" / («^/^) , which, for a 
purely radiation dominated universe (i.e. with 62 = 
and hi 7^ 0), is identically vanishing, whereas the po- 



tential for u is 



{aVT[ 



I 



which in the 



same situation would be ^ 2/r/'^. During the transition 
however, the presence of an amount of matter, however 
tiny, leads to a nonvanishing 62, and hence a nonzero di- 
verging term for small conformal times ~ h\j {^h\r\cqr\): 
the radiation dominated universe represents a singular 
limiting case. This means that both potentials are large 
already at small times, deep into the radiation era, and 
the approximation fc^ <C Vr_m is accurate both before 
and after the transition and for both variables. This ac- 
counts for the fact that the Bardeen potential changes 
during the transition, but only insofar as the amplitude 
is concerned, leaving its spectrum unaltered. This is be- 
cause, in this situation, there is no potential crossing: the 
modes are always below the barrier. Fig.|Hlillustrates this 
fact and summarizes the situation by showing a sketch of 
the perturbation potential 1^_m together with the evo- 
lution of the gravitational potential. 

Another situation of cosmological interest to compare 
the bounce with is a phase of quasi-exponential inflation 
followed by preheating and the subsequent epoch of radi- 
ation domination. When only one field is present, the po- 
tentials for either u or v are essentially undistinguishable 
and both coincide numerically with the inverse Hubble 
size , as shown schematically in Fig. |51 For more than 
one field, the situation is qualitatively different and can- 
not be understood by means of a simple potential ,2^ . 
For a given wave number fc, the spectrum is frozen when 
the wavelength hits the potential, which is often phrased, 
because of the similarity with the Hubble scale, as "hori- 
zon exit" (see Ref. 1151 for a more detailed discussion of 
this point, and Ref. jig in the bounce context). As illus- 
trated in Fig. O the crucial difference between the two 
situations, namely preheating and bounce transitions, is 
that in the latter case the potential and the Hubble scale 
behave in completely different ways whereas they corre- 
spond in the former, at least in the region of potential 
crossing. Far from the bounce itself, however, the poten- 
tial tends to li? again, in a fashion similar to what hap- 
pens during inflation. We conclude that in the bounce 
case, the potential is the quantity that matters and the 
Hubble scale is irrelevant for the calculation of the am- 
plification of perturbations. As a consequence, for prac- 
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FIG. 8: The effective potential Vei-m(?7) for the perturbation 
variables u{Tj) or v{rj) for the radiation to matter transition, 
derived from the scale factor given by Eq. This log-log 

sketch shows the potential (full line) for either of the vari- 
ables (they differ by numerical factors) as well as the exact 
solution (for fe = 0) for the Bardeen potential (dashed line) 
as a function of conformal time. 



tical calculatory purposes, the phrase "Hubble crossing" 
appears misleading in this context p6i | and the phrase 
"potential crossing" should be used instead. 



VI. CONCLUSIONS 

In this section, we summarize the main results ob- 
tained above and discuss them in a more general frame- 
work. 

Assuming general relativity as the theory describing 
gravitation during a bouncing stage happening in the 
early universe, letting the matter content be in the form 
of a scalar field, and restricting attention to the closed 
spatial section case in order to satisfy the null energy con- 
dition, we were able to develop a general formalism by 
expanding any bouncing scale factor around the /C = 1 de 
Sitter-like bouncing solution. This expansion is charac- 
terized by two parameters S and ^ which, in some sense, 
are the counterparts of the slow-roll parameters in the 
usual inflationary models '24]. Because this expansion 
permits a general calculation of the potential for the pri- 
mordial scalar gravitational perturbations, this allows to 
fully determine the structure of their evolution as they 
propagate across the bounce. 

The potential Vu obtained is radically different from 
the Hubble scale at the relevant times. This has to be 
contrasted with the inflationary paradigm for which TC^ 
and Vu are almost identical. 

An important conclusion of this work is that a bounce 
phase, even a short one, can affect large scales of per- 
turbations. General arguments aiming at showing the 
contrary therefore suffer from our counter-example. The 



bounce itself is part of the mechanism described in the 
Introduction, so that the transfer matrix we obtained 
participates to the one of Eq. ^ through 

lim T cx T.< • k-^ ■ T.^ , (88) 

where the k dependence stems from the solution H78|) and 
the unknown matrices and refer to the unknown 
parts sketched in Fig. ^| The coefficients one is inter- 
ested in, namely Tn and T12, giving the amplitude of the 
growing mode in the expanding phase as functions of the 
modes in the contracting phase, accordingly can depend 
on A;. In addition, it is important to notice that, as shown 
m Ref. Pi, this mechanism does not violate causality; a 
similar statement was also emphasized in Ref. jSflJ. 

Paradoxically, obtaining a spectral modification at the 
bounce is possible provided the bounce lasts the mini- 
mal amount of conformal time compatible with the NEC 
preservation. Nevertheless, the assumption of no effect 
can be justified provided the constraint 7/0 ~ 1 5^ 1 is sat- 
isfied, or in the pure de Sitter case having rjo — I strictly. 
This last situation is what happens in models in which 
the bounce takes place for a vanishing value of the scalar 
field kinetic energy , whereas the former case implies a 
kinetic energy density (not the scalar field itself) for the 
scalar field comparable to the Planck scale, which may 
render the semi-classical field theory dubious. 

This can be particularly important in view of the string 
motivated potential alternatives to inflation of the pre big 
bang kind if it turns out that these models might lead to 
such spectral corrections as discussed above. This condi- 
tion needs be verified in each particular situation. For in- 
stance, in the pre big bang case, one would need to model 
the bounce occurring in the Einstein frame, in which our 
formalism is well suited, to see what the behaviour of 
Vu is in this context. Therefore, and unfortunately, one 
consequence of the failure of any general argument pre- 
venting any alteration of the spectrum is that one needs 
to explicitly model a regime in which higher order string 
corrections are dominant. Avoiding this was the main 
interest of the general argument in question. 

We also obtained that the relevant propagation vari- 
able is not V, whose fiat space equivalent is commonly 
used for quantization, i.e. for setting up the initial con- 
ditions, but rather the intermediate variable u, directly 
related to the Bardeen potential. This is to be compared 
with what was recently obtained in Ref. 0], based on a 
completely different theory of gravity, in which neither 
variable happens to be bounded at the bounce. 

The spectrum of gravitational wave cannot be affected 
by propagating through these bounces. This exemplifies 
the fact that there is no fundamental reason according to 
which scalar and tensor modes should propagate similarly 
through a bounce. 

The picture that emerges for the construction of a com- 
plete model of the universe is shown in Fig. 1101 and con- 
sists in a regime in which quantum field theory in a time- 
dependent background is well suited, as is the case for in- 
stance in many string motivated scenarios 7, 8] ; this first 
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9: Left panel: Effective potential U and inverse horizon size relative to the scale of the perturbations in inflation 



models as functions of the conformal time r\. The inflation phase, in this sketch, is smoothly linked with the radiation dominated 
epoch (RDE). The times at which the effect of the potential is comparable with the scale, i.e. ~ JJ are seen to be essentially 
the times at which the scale enters and exits the horizon, i.e. k ^ TL, and are hence labeled "h.c," standing for horizon crossing. 
The primordial power spectrum (PPS) is understood to be the spectrum that is obtained in the phase for which the modes 
are frozen and indicated by an arrow. The actual power spectrum, in such a model, also needs to pass the radiation to matter 
domination transition later on. Right panel: Effective potential Vu and inverse horizon size Ti.^ relative to the scale of the 
perturbations in the bounce model as functions of the conformal time rj. The difference with the inflation case is striking. 




bation that would be sort of pre-primordial. Then, unless 
the curvature was always important in this first period, 
it is followed by an unknown epoch which connects to 
the bounce itself, which should also be followed by yet 
another unknown epoch in order for the curvature to be 
negligible jl6l |. This reveals the most important differ- 
ence between bouncing scenarios and inflation, namely 
the need for a high curvature phase, which we have seen 
may drastically modify the physical predictions. 



FIG. 10: The effective potential Vu{ri) for the perturbation 
variables ii(r;) for our bounce model when one connects this 
bounce transition to both a previous contracting phase on one 
side and to the usual radiation dominated phase later on the 
other side. 



phase allows an easy calculation of a spectrum of pertur- 
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